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ABSTRACT 

We present HST STIS observations of the galaxy NGC 4382 (M85) and axisymmetric models of the galaxy 
to determine mass-to-light ratio (Tv) and central black hole mass (Mbh)- We find Ty — 3.74 ±0.1 M^.jLp, 
and Mbh = l -3^i 2 x 10^ at an assumed distance of 17.9 Mpc, consistent with no black hole. The upper 
hmit, Mbh < 9.6 x 10^ M0(2c7) or Mbh < 1-4 x 10^(3(j) is consistent with the current M-a relation, which 
predicts Mbh = 8.8 x 10^ M,., at Oe — 182 km s^', but low for the current M-L relation, which predicts Mbh = 
7.8 X 10^ M,, at Lv = 8.9 x 10^" U,y. HST images show the nucleus to be double, suggesting the presence of 
a nuclear eccentric stellar disk, in analogy to the Tremaine disk in M3 1 . This conclusion is supported by the 
HST velocity dispersion profile. Despite the presence of this non-axisymmetric feature and evidence of a recent 
merger, we conclude that the reliability of our black hole mass determination is not hindered. The inferred low 
black hole mass may explain the lack of nuclear activity. 

Subject headings: galaxies: individual (NGC 4382, M85) — galaxies: kinematics and dynamics — black hole 
physics — galaxies: nuclei 



1. INTRODUCTION 

Finding a black hole at the center of a galaxy is no longer 
a surprise. The prevalence of these black holes is well es- 
tablished (Richstone et al. 1998). Their importance has also 
been recognized, for example, as active galactic nuclei cen- 
tral engines (Rees 1984). The tight correlation of black hole 
masses with host galaxy properties strongly suggests an un- 
derlying link between galaxy and black hole evolution. The 
black hole mass has been found to be correlated with the stel- 
lar spheroid's mass (Dressier 1989; Magorrian et al. 1998; 
Laor 2001; McLure & Dunlop 2002; Marconi & Hunt 2003; 
Haring & Rix 2004), luminosity (the M-L relation, Kormendy 

' Based on observations made with the Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with GO pro- 
posals 5999, 6587, 6633, 7468, and 9107. 



1993; Kormendy & Richstone 1995; Magorrian et al. 1998; 
Kormendy & Gebhardt 2001; Gultekin et al. 2009c), stellar 
velocity dispersion (the M-a relation, Ferrarese & Merritt 
2000; Gebhardt et al. 2000; Tremaine et al. 2002; Gultekin 
et al. 2009c), galaxy core parameters (Milosavljevic et al. 
2002; Ravindranath et al. 2002; Graham 2004; Ferrarese et al. 
2006; Merritt 2006; Lauer et al. 2007b; Kormendy & Ben- 
der 2009), and globular cluster system (Burkert & Tremaine 
2010). Current theoretical work addressing these scaling re- 
lations focuses on feedback from outflows that are powered 
by accretion onto the central black hole. Thus the conse- 
quences of black hole growth plays a role in regulating star 
formation in the galaxy (e.g.. Bower et al. 2006; Hopkins et al. 
2006). Despite this progress, we do not have a complete un- 
derstanding of the physics involved, and the pursuit of black 
hole masses remains important. 

Masses found from primary, direct, dynamical measure- 
ments are the basis from which all other black hole masses are 
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derived. Indirect mass indicators, such as AGN line widths, 
are calibrated to reverberation mapping, direct yet secondary 
measurements (Bentz et al. 2006), which are themselves nor- 
malized against the direct dynamical measurements (Onken 
et al. 2004; Woo et al. 2010). Currently, even the empirical 
scaling relations are incomplete. The M-L and M-a rela- 
tions make different predictions at the upper end (Lauer et al. 
2007a). The possibility of increased intrinsic scatter at the 
low end remains untested (Volonteri 2007), and there is grow- 
ing evidence that late-type galaxies as a whole and pseudo- 
bulges in particular do not lie on the same relation as early 
type galaxies (Hu 2008; Giiltekin et al. 2009c; Greene et al. 
2010; Kormendy et al. 2011). There is also much interest in 
multi-parameter extensions to these relations (AUer & Rich- 
stone 2007; Hopkins et al. 2007). 

It is this last open question, the existence and utility of 
multi-parameter scaling relations, for which the galaxy in this 
study, NGC 4382 (M85), was selected. NGC 4382 Hes in a 
narrow range in velocity dispersion (180 < C7 < 220 km s^', 
for which Mbh 10^ M ) based on HyperLEDA^ central ve- 
locity dispersion measures (Paturel et al. 2003). With enough 
galaxies from a narrow range in velocity dispersion, we may 
test for additional trends in black hole mass with other host 
galaxy parameters. This range was chosen because galaxies in 
this range may have either core or power-law surface bright- 
ness profiles and because both late-type and early-type galax- 
ies lie in this range. The galaxies were also selected based on 
their distances so that the predicted radius of influence was 
larger than 0"1. The radius of influence is defined as 



^infl = 



GMbh 

72 (/?,nfl) ' 



(1) 



where the velocity dispersion (y{R) is a function of projected 
distance from the center along the major axis and is evaluated 
at the radius of influence. Black hole masses were estimated 
from their central velocity dispersion measurement and the 
M-a fit due to Tremaine et al. (2002). 

NGC 4382 is an E2 galaxy (Kormendy et al. 2009) with dif- 
fuse stellar light surrounding it, which has led some to classify 
it as an SO (de Vaucouleurs et al. 1991). We take the distance 
to NGC 4382 to be 17.9 Mpc (calculated assuming a Hubble 
constant of = 70 km s^^ Mpc^'). The surface bright- 
ness profile as a function of radius reveals a core at the center 
(Lauer et al. 2005) and may be parameterized with a "Nuker 
Law" given by 
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which is a broken power-law profile with variable sharpness 
in the break (Lauer et al. 1995). NGC 4382 has ry = 0"93 = 
80.7 PC 4 = 15.67 mag arcsec-2, a = 1.13, j3 = 1.39, and 
7 = 0.00 (Lauer et al. 2005). The total luminosity of the 
galaxy is My = -22.54. 

In section 2 we describe the observations and data reduc- 
tion, including new space spectroscopic observations (section 
2.1), ground-based spectra (section 2.2), and imaging data 
(section 2.3). The kinematic modeling and its results are pre- 
sented in section 3, and we discuss the caveats for and impli- 
cations of our results in section 4. 



2. OBSERVATIONS 
2.1. STIS Observations and Data Reduction 

Measuring black hole masses precisely requires spectra 
with high spatial resolution. Thus, most precise black hole 
mass measurements come from observations using the Hub- 
ble Space Telescope (HST), though adaptive optics techniques 
are showing promise (e.g., Krajnovic et al. 2009; Nowak 
et al. 2010, Gebhardt et al. in preparation) and mass measure- 
ments using maser observations are ramping up (Kuo et al. 
201 1). We observed Ca II triplet absorption from NGC 4382 
with the Space Telescope Imaging Spectrograph (STIS) on 
HST set with the G750M grating and a 52" x 0"2 slit. The 
medium-dispersion grating, as opposed to the low-dispersion 
G750L grating, is necessary in order to get the spectral reso- 
lution high enough to recover line-of-sight velocity distribu- 
tions (LOSVDs) with sufficient precision. The slit was set at 
a width of 0"2 to optimize the signal-to-noise ratio as core 
ellipticals have relatively low central surface brightness and 
widening the slit as far as possible allows as short an obser- 
vation as possible. The slit was positioned along at a position 
angle of PA ~ 48° east of north close to the the photomet- 
ric major axis position angle of PA — 30° as determined from 
//5'7'AVFPC2 observations (see section 2.3 below, Lauer et al. 
2005; Kormendy et al. 2009). We obtained 16 exposures at 
5 dither positions for a total of exposure of 18 794 s. The 
STIS CCD has a 1024x 1024 pixel format, a readout noise of 
^ le^ pix~\ and a gain of 1.0 without on-chip binning. The 
spectra spanned a wavelength range of 8257-8847 A, and our 
wavelength solutions revealed a reciprocal dispersion of 0.554 
A pix"^ and the spatial scale was 0"05071 pix"' for G750M 
at 8561 A. 

The STIS data reduction was done with routines developed 
for this purpose (Pence 1998; Pinkney et al. 2003) follow- 
ing the standard pipeline: Raw spectra were extracted from 
the multi-dimensional FITS file, and then a constant fit to the 
overscan region was subtracted to remove the bias level. The 
STIS CCD has warm and hot pixels that change on timescales 
of about a day. These pixels require that the subtraction of 
dark current be accurate, which we accomplished by using the 
iterative self-dark technique (Pinkney et al. 2003). After flat- 
fielding and dark subtraction, spectra were shifted vertically 
to a common dither, combined, and rotated. One-dimensional 
spectra were then extracted using a bi-weight combination of 
rows. Near the galaxy center, we adopted a 1 -pixel wide bin- 
ning for maximum spatial resolution. These data-reduction 
methods are similar to those in Pinkney et al. (2003), which 
the interested reader may consult for details. 

LOSVDs were extracted from reduced spectra as described 
in Giiltekin et al. (2009b). Each galaxy spectrum is a con- 
volution of the intrinsic spectrum of stars observed in the 
aperture with the LOSVD of those stars. We deconvolved 
the observed galaxy spectrum using the template spectrum 
composed from standard stellar spectra using a maximum 
penalized-likelihood method (Gebhardt et al. 2003; Pinkney 
et al. 2003). We present Gauss-Hermite moments of the ex- 
tracted velocity profiles in Table 1 and Figure 1. Although 
it is common practice to communicate velocity profiles with 
Gauss-Hermite moments, we use LOSVDs binned in velocity 
space for our modeling described below in section 3. 



^ Available at http://leda.univ-lyonl.fr/. 



2.2. Ground-Based Spectra 



BH in NGC 4382? 



3 



Table 1 

Kinematic Profile for NGC 4382 from STIS Observations 



R 


Width 


V 


a 


h 


/14 


(") 


(pix) 


(kms"') 


(kms"') 






0.00 


1 


15±17 


147±20 


-0.017±0.06 - 


-0.068±0.060 


0.05 


1 


4±14 


162±12 


-0.049±0.05 - 


-0.045±0.036 


0.10 


1 


24±15 


150±13 


-0.006±0.05 - 


-0.067±0.029 


0.18 


2 


31±15 


145±12 


-0.071±0.05 - 


-0.049±0.032 


0.30 


3 


14±16 


170±17 


0.021±0.06 - 


-0.063±0.046 


0.58 


8 


24±18 


178±20 


-0.021±0.07 - 


-0.016±0.053 


1.12 


14 


20±19 


166±20 


-0.001±0.07 


0.015±0.060 


Note. 


— Gauss-Hermite moments for velocity profiles 


derived from 



STIS data. First and second moments are given in units of km s" . Radii 
are given in arcsec, and the second column gives the width of the radial bin 
in pixels, which are 0"05 1 . 

Ground-based velocity information was obtained from 
archival data using Spectrographic Areal Unit for Research on 
Optical Nebulae (SAURON)^, an integral-field spectrograph 
unit mounted on the William Herschel Telescope in La Palma 
(Bacon et al. 2001). As the SAURON instrument and observa- 
tions were designed to measure and characterize the internal 
kinematics of the selected galaxies (de Zeeuw et al. 2002), the 
data are excellent for our purposes. In its low resolution mode, 
the instrument has a field of view of 33 x 41" with 0"94 pix- 
els, each of which provides a spectrum with FWHM = 4.2 A 
spectral resolution. This spectral resolution corresponds to 
CTinstr = 108 km s"^ at 4950 A, near the center of the 4800- 
5380 A wavelength range. The data were taken on 14 March 
2001 in two pointings, each consisting of 4 exposures of 1800 
s each under FWHM= 2"7 seeing (Emsellem et al. 2004). 

The stellar kinematics data are provided as Gauss-Hermite 
moments for each lenslet position on the galaxy. The Gauss- 
Hermite moments from each lenslet were converted into 
LOSVDs. The data from each quadrant of the galaxy were 
combined into one, changing signs of the odd moments as 
necessary. We binned the data into 4 position angles (0, 20, 
30, and 70° east of the major axis) with 7 radial bins and one 
position angle (45° east of the major axis) with 6 radial bins 
for a total of 34 LOSVDs from ground-based data. 

To incorporate errors in the Gauss-Hermite moments in our 
data set we created lO'* Monte Carlo realizations of LOSVDs 
for each SAURON lenslet. As described in Gultekin et al. 
(2009b), Gauss-Hermite moments corresponding to unphys- 
ical negative values are assigned a value of zero with a con- 
servative uncertainty. For each spatial bin, we took the me- 
dian LOSVD of all realizations at a given velocity and the 
standard deviation as the error, which dominated the individ- 
ual measurement errors. The Gauss-Hermite moments of the 
LOSVDs are presented in Figure 1 as a function of radius 
along with our STIS data. These measurements agree well 
with the ground-based kinematics from Fisher (1997). 

The profiles were binned into 1 3 equal bins in velocity from 
—500-500 km s^' about the systemic velocity covering the 
range of velocities measured. Using the ground kinematic 
data, we compute an effective stellar velocity dispersion, de- 
fined as 

._ Io' {a{rf + V{rf)l{r)dr 

.k nndr 

where Re is the effective radius, I{r) is the surface brightness 



profile (see section 2.3 below), and V (r) and a{r) are the first 
and second Gauss-Hermite moments of the LOSVD. Using a 
non-parametric method of integrating the brightness and el- 
lipticity profile Kormendy et al. (2009) find = 102 ±6". 
From the ground-based velocity profile, we find an effective 
stellar velocity dispersion of C7^ = 182 ± 5 km s^ ' . 

2.3. Imaging 

The high-resolution photometry of the central regions of 
NGC 4382 comes from Wide Field Planetary Camera 2 
(WFPC2) observations on HST using filters F555W (V) and 
F814W (/). The observations, data reduction, and surface 
brightness profiles (including Nuker profile fits) are detailed 
by Lauer et al. (2005). Surface brightness profiles are also 
available at the Nuker web page."* The wide-field data that 
we use are literature data from ground observations on the 
1 .2-m telescope of the Observatoire de Haute-Provence with 
PSF FWHM== 3" 12 originally obtained by Michard & Mai-- 
chal (1994). The ground data are B-band data, which we have 
color-corrected with B — V = 0.9, determined by requiring the 
space data and ground data to match where they overlap in the 
spatial direction. Figure 2 shows the surface brightness pro- 
files of the space and ground data along with our adopted, 
combined surface brightness profile for our modeling. 

Since our modeling efforts began, a superior ground-based 
photometric data set became available (Kormendy et al. 
2009). The data come from several different instruments. In- 
side of 1", it is the same WFPC2 F555W data that we use. 
Starting at 1", the profile is an average of WFPC2 and ACS 
data as well as two different cameras on the Canada-France- 
Hawaii Telescope. At the largest radii, the data come from the 
McDonald 0.8-m telescope, which has a wide field of view. 
We plot it in Fig. 2 for comparison. The agreement inside of 
R < 15" is striking. At large radii, there is a small but sys- 
tematic deviation from our adopted surface brightness profile. 
This is almost certainly due to the different colors used and 
is likely evidence of a slight color gradient starting outside of 
R > 100". At these large radii, the smaller PSF of the Kor- 
mendy et al. (2009) data set does not have a large advantage, 
and the extremely wide coverage — out to R = 625" — is not 
used since our kinematic profiles only go out to R = 30". 

3. KINEMATIC MODELING 

We use the three-integral, axisymmetric Schwarzschild 
method to make kinematic models of NGC 4382. The model 
is constructed in several steps. First the observed surface- 
brightness profile, E(r), is deprojected into an axisymmetric 
luminosity density, j{r,9). This deprojection assumes that 
surfaces of constant luminosity density are coaxial ellipsoids 
and depends on a chosen inclination, which we take to be 
/ = 90°. Under the assumption of constant mass-to-light ra- 
tio of a chosen value, T, the mass density of stars is trivially 
obtained p (r, 0) = T j. We can then calculate the stellar gravi- 
tational potential from Poisson's equation. The potential from 
a point mass with a chosen mass, Mbh, is then added at r = 0. 

With the potential for the entire system in hand, we then 
calculate orbits of representative stars. The number of orbits 
is different for each mass model with more orbits needed to 
sample the phase space when T or Mbh is large, but the num- 
ber ranged from 15735 to 17213 orbits. The number of orbits 
is increased by increasing the density of orbits in energy, an- 
gular momentum, and non-classical third integral space. We 



^ Downloaded from http://www.strw.leidenuniv.nl/sauron/ in Feb 2010. 



^ See http://www.noao.edu/noao/staff/lauer/wfpc2_profs/. 
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Figure 1. Gauss-Hermite moments of LOSVDs for NGC 4382 as a function of radius, from top to bottom: V, fj, /13, and h4. Values at tlie center of tlie 
galaxy (R = 0) are plotted at R = ()!'025. Moments from opposite sides of the galaxy are symmetrized and combined. The sense of the sign is such that radius 
is increasing to the southwest. Crosses are Gauss-Hermite moments of LOSVDs from HST STIS data along the major axis. Also plotted are ground-based 
Gauss-Hermite moments of LOSVDs along the major axis (black diamonds) and skew axes (other colors with PA relative to the major axis indicated on the right 
in the top panel) from SAURON archives. Though Gauss-Hermite moments are not fit directly in the modeling, the jagged lines are the resulting Gauss-Hermite 
fit to the best-fit model's LOSVDs from section 3 for the major axis (solid) and skew axes (dashed). The best-fit model has Mbh = 1.5 x lO' M, and Yy = 3.72. 
We indicate the radius of the sphere of influence of a black hole with mass predicted by the M-a relation (light gray) or with our best-fit mass (medium gray). 
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Figure 2. Surface brightness of NGC 4382. Red and green crosses (upper 
and lower) show F814W and F555W WFPC2 data, respectively, from Lauer 
et al. (2005). Blue circles show Michard & Marchal (1994) ground photom- 
etry from the 1.2-m telescope of the Observatoire de Haute-Provence with a 
PSF FWHM of 3!' 12. Green circles are the same data but color-shifted by 
0.9 magnitudes to match the F555W data. The solid curve is our assumed 
V-band surface brightness profile of NGC 4382. The purple line is data from 
Kormendy et al. (2009) and comes from a variety of instruments and bands. 
The difference between our adopted surface brightness profile and the puiple 
line is small and can be attributed primarily to slightly different colors at large 
radii. 



increase the number of orbits with sparsely gridded parame- 
ter space (Menand Ty) until our results converge, and then 
we run a finer grid in parameter space, which we report here. 
The amount of time each representative star spends in a given 
bin in position and velocity space is monitored so as to iden- 
tify the orbits' possible contribution to the observed surface 
brightness and LOSVD. For each set of parameters, we find 
the non-negative weights and a goodness-of-fit (j^) for the set 
of orbits that best matches the observed LOSVD while repro- 
ducing the observed surface brightness. The entire method is 
explained in more detail in Gebhardt et al. (2003) and Siopis 
et al. (2009). 

To determine what values of Mbh and Tv to consider, we 
ran a sparse grid of a wide range of values and then re- 
fined around the best-fitting initial guesses. The final ranges 
considered were Ty — 3.5-4.0 M,/L;, v and Mbh = to 
2x10^ The sampling of parameters was not uniform but 
may be ascertained from the dots in Figure 3. This figure also 
plots contours of in the Mbh-^v plane. The best-fit model 
has parameters Mbh = 1 -5 x 10^ Mq and Ty = 3.72 M^jL^y, 
but models with Mbh = are consistent at about the 1 a level. 
The Gauss-Hermite moments of the LOSVD for the best- 
fit model are shown in Figure 1 and show good agreement 
with the data. To obtain our final estimates for Mbh and 
Tv, we marginalize over the other parameter to get Mbh = 



1.3 



+5.2 
1.2 



10^ Mj, and = 3.74 ±0.10 M; /i,y, where the 



errors are Iff uncertainties derived from changes from the 
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Black Hole Mass (10* M,) 

Figure 3. Mass modeling contours for NGC 4382. Contours are for 
= 1.00, 2.71, 4.00, and 6.63, which bracket individual parameter con- 
fidence levels of 68.3, 90.0, 95.4, and 99.0%, respectively. Contours have 
been smoothed for plotting, and each level is filled with a solid color The 
square shows the best-fit model. Dots indicate parameters modeled. The 
best-fit model has Mbh = 1.5 x lO'' Mg and Tv = 3.72 Mg/Le.v- Marginal- 
izing over the other parameter we find Mbh = 2 ^ ™d = 
3.74±0.10A/;/L,;,,v. 

minimum marginalized of 1^%^ = 1 ■ Thus, we cannot rule 
out the absence of a black hole, but we may put 2 and 3 a 
upper limits at Mbh < 9.6 x 10^ and Mbh < 1-4 x 10^ M?, 
respectively. We plot marginalize Aj^ as a function of Mbh 
and Tv in Figure 4. 

Figure 5 shows the velocity dispersion tensor for the best- 
fit model by plotting the ratio of the radial velocity disper- 
sion (Or) to the tangential velocity dispersion (o;), defined as 
of = {Oq + o'^)/2, where is the second moment of the az- 
imuthal velocity relative to the systemic velocity rather than 
relative to the mean rotational speed. An isotropic velocity 
dispersion would have C7,./(7, = 1. Uncertainties are ^ 0.2, 
estimated from the smoothness of the profiles (Gebhardt et al. 
2003). There is a clear transition at 0"1. Outside of this ra- 
dius, the orbits are roughly isotropic (dr/Of ~ 1), but inside 
of this radius the orbits are tangentially biased (cT,-/<7r < 1). 
There is a very strong indication of rotation at the center of 
the model galaxy along the major axis ((7^/ Oi — 0.24). 

4. DISCUSSION 

Before discussing the implications of our results, several 
potential difficulties in interpretation should be mentioned for 
this galaxy. Our kinematic modeling requires three assump- 
tions: (1) a constant stellar mass-to-light ratio over the range 
of interest of the galaxy with the exception of the central black 
hole, (2) axisymmetric mass distribution, and (3) that the sys- 
tem is in dynamical equilibrium. We discuss each of these in 
turn. 

4.1. Constant T? 

The first of these assumptions, constant mass-to-light ra- 
tio, is unlikely to be far from reality. Using XMM-Newton 
and Chandra data, Nagino & Matsushita (2009) study the 
gravitational potential as revealed by X-ray emission from 
the interstellar medium, assumed to be in hydrostatic equi- 
librium and spherically symmetric. Within ^ 2 kpc, about the 
outer extent of our data, they find a constant B-band mass- 
to-light ratio consistent with a potential dominated by stellar 




Mass-to-light ratio 



Figure 4. Marginalized H^X' curves as a function of Mbh (top) and Tv (bot- 
tom). Dashed lines show intersections of the curves with = 1,4, and in 
the case of the top panel 9, which correspond to the 1, 2, and 3o' confidence 
levels. There is a slight non-monotonicity in the top panel as one increases 
in Mbh away from the minimum, but it is small and entirely within the 1 a 
interval. Because the deviations from monotonicity are small (i.e., less than 
unity), we may be confident that they are not significant. In both cases there 
is a clear minimum near the best-fit value of each parameter, but in the case 
of Mbh, the data are consistent with Mbh = at about the IcT level. 



mass. Over the range of 0"2 <r < 10" (« 0.2-1 kpc) NGC 
4382 has a shallow V — I color gradient d{V — I)/d\og{r) — 
+0.006 ± 0.003, i.e., growing bluer with decreasing radius, 
with y - / = 1.108 ±0.001 at r = 1" (Lauer et al. 2005). 
While the sign of the color gradient is unusual, the magni- 
tude is small enough to dismiss worries about changing stellar 
population in the galaxy. 

To quantify the effect of deviations from constant mass-to- 
light ratio on estimates of black hole mass, we note that when 
the black hole kinematic influence is barely resolved, Mbh 
and T are anti-correlated (Gebhardt & Thomas 2009; Schulze 
& Gebhardt 2011). Thus if there were unaccounted system- 
atic effects such as a radial gradient in T, it would increase the 
uncertainty in Mbh- We can easily estimate the magnitude of 
the effect since it is linear with T. That is, a 30% change in T 
will result at most in a 30% change in Mbh inference. Based 
on the radial color gradient analysis above, the magnitude of 
any radial T gradient must be small, a couple of percent at 
most. Thus the total systematic uncertainty in Mbh is less 
than a couple of percent. 
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Figure 5. Shape of the velocity dispersion tensor for NGC 4382 from the 
best-fit model orbit-solution. The black line is along the major axis, and the 
other lines show skew angles with position angles relative to the major axis 
as given in the legend. The values for the central part of the galaxy are plotted 
at a radius of 0!'01. The dashed line shows the radial extent of the ground- 
based spectroscopic data. The orbits are isotropic until the very center at 
which point they become tangentially biased. We indicate the radius of the 
sphere of influence of a black hole with mass predicted by the M-a relation 
(light gray) or with our best-fit mass (medium gray). 

4.2. Axisymmetric or eccentric disk? 

The second of these assumptions, axisymmetric mass distri- 
bution, is potentially violated. Note that the counter-rotating 
"kinematically decoupled core" (KDC) found in the OASIS 
data does not necessarily imply deviations from axisymme- 
try (McDermid et al. 2004). As can be seen from the im- 
age of the galaxy (Fig. 7), the change in ellipticity near the 
center (e — 0.6 at r w 0"2) to the outer regions (e = 0.2 for 
r > 1") may demonstrate that this system is, overall, triaxial 
(Chakrabarty 2010). Modeling a truly triaxial system requires 
a triaxial code (van den Bosch & de Zeeuw 2010), but the case 
for triaxiality is not straightforward because inside of 0"2, the 
ellipticity decreases again (Lauer et al. 2005). 

An alternative interpretation of the two-dimensional pho- 
tometry is that there are two peaks in surface brightness at the 
center of the galaxy with projected separation 0"25 (Lauer 
et al. 2005). Such a double nucleus appears to be similar to the 
double nucleus in M31 (NGC0224; Lauer et al. 1993, which 
Tremaine (1995) explained as the result of a projection of a 
central eccentric disk of stars that is stable only in the pres- 
ence of a massive dark object, such as a black hole. 

This interpretation is born out in the unsymmetrized STIS 
spectroscopy (Fig. 8). There is a prominent increase in the ve- 
locity dispersion at the location of secondary surface bright- 
ness peak R k, 0"4. This would be expected for an eccentric 
stellar disk. 

An eccentric disk cannot be modeled by our orbit superpo- 
sition code, which forces axisymmetry. In order to estimate 
the error introduced by forcing axisymmetry, we compared 
the mass that would be obtained by treating an eccentric disk 
as though it were a circular orbit. We use the mass estima- 
tor {v\y) averaged over the orbit, where \\ is the line-of-sight 
velocity and y is the distance from the center of light of the 
orbit on the sky. This is a suitable surrogate for our modeling. 
We then compare the value of the estimated mass assuming 
an elliptical orbit to the mass inferred when assuming a circu- 
lar orbit with the same semi-major axis. Treating an elUptical 




Figure 6. Position of STIS slit from the acquisition camera. This image 
comes from the STIS acquisition image taken just before the slit and dis- 
persive element were added. We have subtracted a symmetric model for 
NGC4382 in order to show the non-symmetric features. The red box shows 
the location of the slit, and the red cross shows the center position of the 
galaxy. Up is to the southwest and positive R in Figure 8. The figure shows 
that the slit does lie on the surface brightness peak. 



orbit as though it were a circular orbit produces an error in the 
black hole mass that depends on the orbit's eccentricity and 
orientation. For an eccentricity of 0.60, the error in the black 
hole mass estimate caused by ignoring the eccentricity varies 
from — 15% to +7%. Orbits with smaller eccentricity produce 
smaller errors. Given the statistical uncertainties in our Mbh 
estimate, we can safely ignore this systematic uncertainty. 

A consequence of our assumption of axisymmetry is that 
we are only sensitive to the presence of a black hole at the 
center of the galaxy. There are two potential issues. The first 
is whether an existing black hole is located at the center; the 
second is whether the STIS slit was positioned to allow mea- 
surement of any central black hole. We consider both of these. 

Because NGC 4382 is a recent merger (see section 4.3), 
a relevant question is whether any black hole in the galaxy 
is located at the center, a requirement for detection with our 
method. The timescale for a black hole to sink to the center 
of a galaxy depends on where it starts its decent and in what 
orbital configuration. In a purely tangential orbit at ^ 20 kpc, 
it would take longer than a Hubble time to rest at the center, 
but a radial orbit at ^ 100 pc would take less than ^ 50 Myr. 
In the most likely case, before the merger, there would be a 
black hole in the primary galaxy. Following the merger, it is 
unlikely that it would have traveled farther than ^100 pc, and 
we can expect any such black hole to be at the center now. 

Since the center of NGC 4382 is morphologically complex 
and we have not found strong evidence for a black hole at its 
center, it is important to determine that the slit was positioned 
in such a way that the central kinematics could be determined. 
From the STIS acquisition camera image (Fig. 6), it can be 
clearly ascertained that it does so. This image comes from 
the STIS acquisition image taken just before the slit and dis- 
persive element were added. We have subtracted a symmetric 
model for NGC4382 in order to show the non-symmetric fea- 
tures. The red box shows the location of the slit, and the red 
cross shows the center position of the galaxy. 

4.3. Dynamical equilibrium? 

The third of these assumptions, dynamical equilibrium, is 
the most important. Our analysis, like any dynamical analysis 
of this kind, does not handle large deviations from dynamical 





Figure 7. Images of NGC 4382. The left panel shows the Kormendy et al. (2009) wide-field image from SDSS gri with a high-pass filter applied to bring out the 
fine structure, oriented such that north is up and east is left. The center panel is from F555W HST/WFPC image and shows the position of the slit with regards to 
the major axis. The right panel is the Lauer et al. (2005) contour of the same F555W ffSTAVFPC data zoomed in to show the double-nucleus. The fine structure 
and double nucleus are evidence of a possible recent merger in this system. The double nucleus is likely a stellar eccentric nuclear disk in analogy to M3 1 . 
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Figure 8. Velocity and velocity dispersion measurements from STIS from 
both sides of the center. The black line shows the best-fit axisymmetric 
model. There is a noticeable increase in C7 at R = +0"5. which corresponds 
to the location of the secondary suii'ace brightness peak, indicative of an ec- 
centric stellar disk or torus. We indicate the radius of the sphere of influence 
of a black hole with mass predicted by the M-fJ relation (light gray) or with 
our best-fit mass (medium gray). 



equilibrium. Schweizer & Seitzer (1992) found a high fine 
structure index in NGC 4382 of E = 6.85, defined as 



I = 5 + log(l+n)+y + B+X, 



(4) 



where 5 is a visual estimate of the strength of the most promi- 
nent ripples with range S — 0-3, n is the number of detected 
ripples, J is the number of jets, B is a visual estimate of the 
maximum boxiness of isophotes, with range B = 0-3, and 
X = or 1 indicates the absence or presence of an X-structure, 
respectively (Schweizer et al. 1990). 

Lauer et al. (2005) consider NGC 4382 to be an excellent 
candidate for a recent merger based on its high value of E, 
one the three-highest non-merging in the Schweizer & Seitzer 
(1992) sample, the KDC mentioned above, and the galaxy's 
blue average V — / color. Kormendy et al. (2009) note that the 
strong fine-structure features are also seen in the gri SDSS im- 
age and interpret this as evidence that the galaxy has recently 
merged but has not fully relaxed. 

The fine structure index, however, is sensitive to gas-rich 
minor mergers that are unlikely to disturb the entire galaxy 
or the central parts of the galaxy from an equilibrium condi- 
tion. Additionally, fine structure may be a poor measure of 
recent gas-poor major mergers (van Dokkum 2005; Tal et al. 
2009). In the case of NGC 4382, the fine structure appears to 
be "Malin shells," which can be the result of gas-poor, primar- 
ily stellar, mergers (Malin & Carter 1980; Quinn 1984). As a 
core elliptical this galaxy is expected to be a mostly gas-poor 
merger (Lauer et al. 2007b; Kormendy et al. 2009). 

The virialization time at the center of the galaxy is much 
shorter than at the outer regions where the Malin shells are 
at apopause and may persist for ^10^ yr so that virialization 
proceeds from the inside out. It is thus not likely to affect the 
inner portions of the galaxy important for our study (Navan^o 
1990). As a final comment on evidence of recent merger ac- 
tivity, we note that Chung et al. (2009) investigated NGC 4382 
for signs of interaction motivated by the ^ 35 kpc projected 
distance from NGC 4395. They found no signatures of inter- 
action in H I or other wavelengths. 

So we urge caution in over-interpretation of our kinematics 
modeling, but for the remainder of this discussion, we take 
our results at face value while considering their implications. 

4.4. Resolution of the sphere of influence 

Resolving the sphere of influence of the black hole is desire- 
able but not necessary for measuring the mass (Giiltekin et al. 
2009c). As the sphere of influence is more and more poorly 
resolved, the uncertainties in mass increase until the estimate 
of black hole mass is consistent with zero or even negative 
(e.g., NGC 3945 Gultekin et al. 2009b). Using our mea- 
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sured effective velocity dispersion of — 182 km s , the 
black hole mass predicted from the M-a relation is Mbh = 
8.8 X 10^ M. . From equation (1) the predicted sphere of in- 
fluence for this object at a distance of 17.9 Mpc has radius 
0"13 or diameter Q'.'26. The central resolution element of our 
STIS spectroscopy is Q'.'l across the slit and 0"05 along the 
slit. Thus the projection of the predicted sphere of influence 
fits entirely inside of the central resolution element. Based on 
our best-fit mass estimate of Mbh = 1.3 x 10^ M,.,, the radius 
of the sphere of influence is 0"03 (diameter 0"06) where the 
the velocity dispersion is 145 km s^^ The diameter of the 
actual sphere of influence is thus about the size of the pixel 
but 3.3 times smaller than the slit width. At such low resolu- 
tion, it is typical to find an upper limit to the mass of the black 
hole, as we have here. The la upper bound to the black hole 
mass is Mbh — 6.5 x 10^ M,,, which has sphere of influence 
radius 0"13, the same as the predicted sphere of influence. 

4.5. Is the black hole mass anomalously low? 

Our modeling reveals that the kinematic observations are 
consistent with no black hole in NGC 4382 at about the la 
level. If there were, in fact, no black hole in NGC 4382, it 
would be one of only two early-type galaxies with upper lim- 
its on their black hole mass below that predicted by the scal- 
ing relations. The other is NGC 3945, but the double bars in 
that galaxy may prevent reliable black hole mass estimation 
(Giiltekin et al. 2009b). Given that our mass measurement 
does not completely rule out the presence of a black hole, we 
consider whether NGC 4382 is a unique object or merely lies 
along the low -end tail of the distribution of black hole masses 
for a given host galaxy property. 

Based on the M-a relation (Mbh/Mq = 
10^'2(aj200 kms-i)'^-24; Giiltekin et al. 2009c), the 
mean logarithmic black hole mass for a galaxy with 
a = 182 kms^' is log(MBH/M.,) = 7.95, corresponding 
to Mbh = 8.8 x 10^ M, . The scatter in the M-a relation is 
log-normal with standard deviation eo = 0.31 ± 0.06 for the 
population of ellipticals. Taking the scatter into account, the 
68% confidence interval of black hole mass is Mbh = 4. x 10^ 
to 1.8 X 10** M:,, which is consistent within our la mass 
estimate. 

To calculate the mass expected from the M-L relation 
(Mbh/M: = 10** ''^(Lv/10"Lv,,v)' "; Gultekin et al. 2009c), 
we must first find the luminosity of the bulge. The to- 
tal luminosity of the galaxy may be obtained from My = 
-22.54±0.05 (Kormendy et al. 2009) using log(Ly/L:,,v) = 
0.4(4.83 - My) (cf. Verbunt 2008). Given that we are adopt- 
ing the more modern classification of E2, this galaxy is "all 
bulge," and so the luminosity is Ly = 8.9 x lO'" A- y- The 
mean logarithmic mass for such a galaxy is log (Mbh/M,;,) = 
8.89 or Mbh = 7.8 x 10** M.... Taking the Eq 0.38 ±0.09 
scatter in the relation into account, the 68% confidence inter- 
val in mass is Mbh = 3.0 x 10** to 2.1 x lO'^ M, . The lower 
limit of this range is more than a factor of 2 larger than our 3 a 
upper limit for the black hole mass. In this sense, the black 
hole mass is anomalously low. 

In order to be consistent with the M-a reltaion but low 
according to the M-L relation, NGC 4382 cannot lie on the 
mean Faber & Jackson (1976) L-a relation. In fact, core 
galaxies with My = —22.54 have a mean logarithmic velocity 
dispersion corresponding to a = 253 km s^', compared our 
measured value of Oe — 182 km s^'. 



The stellar mass of the bulge is M* — LyTy = 3.3 x 
lO'^ Mg. This is consistent with g- and z-band model mag- 
nitudes, which give M* 4.0 x lO" M; (Befl et al. 2003; 
Gallo et al. 2010). Thus Mbh/M* = 3.9 x 10"^ compai-ed to 
the standai-d 1.3 x 10"^ value (Kormendy & Gebhardt 2001). 
Using the 3a upper mass bound, the ratio is Mbh/M* = 
4.2 X 10^"*, still a factor of 3 lower In section 4.6 we dis- 
cuss how our value of Ty may be wrong because we do not 
include a dark matter halo in our modeling, but this is only 
likely to result in a small (5-10%) decrease in Ty (Gebhardt 
& Thomas 2009; Schulze & Gebhardt 2010). 

Like many other early type galaxies, NGC 4382 has a 
flat luminosity core. It has been argued that the cores are 
scoured out by the inspiral of black holes during galaxy merg- 
ers (Begelman et al. 1980; Ebisuzaki & Makino 1995; Faber 
et al. 1997; Volonteri et al. 2003a,b; Milosavljevic & Merritt 
2003; Lauer et al. 2007b; Kormendy & Bender 2009). This 
process ejects stars on elongated orbits and leads to tangen- 
tially biased stellar distribution functions (See Fig. 5). The 
galaxy's core can be described by its stellar mass "deficit," the 
mass in stars ejected from what was previously a power-law 
profile. The most recent numerical simulations on black hole 
mergers find that for nearly equal masses, the mass deficit 
should scale with the total mass of the binary (Merritt 2006); 
but for mass ratios far from unity, the mass deficit should scale 
with the mass of the secondary (Sesana et al. 2008). 

Kormendy & Bender (2009) included NGC 4382 in their 
study of the correlations between black hole mass and mass 
deficits. They estimated Mbh = 1.0 x 10** M, by assum- 
ing that it follows the Tremaine et al. (2002) M-a relation, 
and they estimated Ty =8.3 by assuming that it follows 
Ty Ly^^ due to Cappellari et al. (2006). Using the photome- 
try from Kormendy et al. (2009), they measured a luminosity 
deficit of 1.6 X 10^ A- y- For their assumed values of Mbh 
and Ty, NGC 4382 lies along the same correlation between 
Mbh and mass deficit as the rest of their sample (roughly 
■'Wdeficit — IOMbh)- Using the values for Mbh and Ty we mea- 
sure here, NGC 4382 appears to have a small black hole mass 
for its core mass deficit, but is consistent at the la level. For 
our best estimate of the mass, Mdeficit/MBH = 45.6, and for the 
la upper bound of the mass Mdeficit/MBH =9.1. Kormendy 
& Bender (2009) also found that the ratio of light deficit to 
total bulge luminosity (Ly.deficit/^'V = 1.8 x 10"^ for NGC 
4382) scales with the ratio of black hole mass to bulge stel- 
lar mass. Using their empirical correlation and our measured 
value of Mbh/M* predicts a ratio of Ly.deficit/i'V = 2.0 x 10^"*. 
With our la upper bound on Mbh, it predicts iy.deficit/i'V = 
1 .7 X lO^-', consistent with the observed value within the total 
scatter of the relation. So it appears that, in sum, the mass of 
any black hole in NGC 4382 is consistent with M-a and core 
scaling properties but low based on M-L. 

The low black hole mass does, however, clear up an appar- 
ent disparity in nuclear activity. The nuclear X-ray luminosity 
between 0.3 and 10 keV is less than Lx < 2.7 x 10^^ erg s"' 
(Sivakoff et al. 2003; Gallo et al. 2010). The core was also 
not detected in radio, typically considered a proxy for jet ac- 
tivity from an accreting black hole. With an X-ray limit and 
a radio detection, a black hole mass limit can be estimated 
(Giiltekin et al. 2009a). Very Large Array observations at 8.4 
GHz, however, did not detect any core radio emission with 
a 3a upper limit of Fcok < 0.11 mJy (Capetti et al. 2009). 
Capetti et al. (2009) note that it is concerning that a galaxy as 
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large as NGC 4382 shows no sign of nuclear activity. Their 
argument is that, under the assumption that NGC 4382 hosts 
a Mbh ^10^ M;, SMBH as would be expected from galaxy 
properties, the absence of nuclear activity in X-ray and radio 
bands indicates a failure of AGN activity to trace black holes. 
Even a small amount of ambient gas could produce enough 
radio emission to be visible for a large black hole. Our black 
hole mass estimate, however, indicates that it is more likely 
that NGC 4382 either has no black hole or a black hole that 
is far undermassive for its host galaxy properties. Therefore, 
the lack of any nuclear activity accurately traces the small or 
non-existent black hole in line with predictions for and exist- 
ing observations of the AGN fraction in Virgo cluster galaxies 
(Volonteri et al. 2008; Gallo et al. 2008). 

4.6. Mass-to-light ratio 

Although this work does not unambiguously detect a black 
hole, the mass-to-light ratio was determined precisely, = 
3.74 ±0.10 M.,/Lr.,y owing to the good coverage and high 
quality of the S AURON data. Unfortunately, NGC 4382 was 
not one of the galaxies modeled by Cappellari et al. (2006); 
so no direct comparison can be made. Note that we do not 
explicitly include a dark matter halo in our modeling. This 
is unlikely to affect our black hole mass results on account 
of the high spatial resolution provided by HST/STIS (Geb- 
hardt et al. 2011). The stellar mass-to-light ratio, however, is 
likely to be affected by not including this component. While 
the inferred stellar mass-to-light ratio changed by a factor of 
2 when including or omitting a dark matter component in the 
modeling for M87 (Gebhardt & Thomas 2009), the high spa- 
tial resolution in our data sets is likely to result in only a 
small (5-10%) decrease in Ty (Schulze & Gebhardt 2010). 
Galaxies with My = —22.54 have a mean mass-to-light ratio 
Yy = 6.7 (Lauer et al. 2007a). Based on the scaling relation 
between Y/ and stellar velocity dispersion in ellipticals (Cap- 
pellari et al. 2006), the implied typical value for an galaxy 
with Ge = 182 km s^^ is roughly Ty w 5.3. The low value of 
Ty in NGC 4382 compared to established ellipticals suggests 
ongoing star formation and would be expected for a recent 
merger. 
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